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Abstract 

We study the evolution of a one dimensional model atom with 5- 
function binding potential, subjected to a dipole radiation field E(t)x 
with E{t) a 27r/a;-periodic real-valued function. Starting with xp{x,t = 0) 
an initially localized state and E{t) a trigonometric polynomial, complete 
ionization occurs; the probability of finding the electron in any fixed region 
goes to zero. 

For %p{x,Q) compactly supported and general periodic fields, we con- 
struct a resonance expansion. Each resonance is given explicitly as a 
Gamow vector, and is 2n/uj periodic in time and behaves like the expo- 
nentially growing Green's function near x = ±00. The remainder is given 
by an asymptotic power series in t~^^'^ with coefficients varying with x. 

1 Introduction 

The ionization of an atom by an electromagnetic field is one of the central 
problems of atomic physics. There exists a variety of approximate methods for 
treating this problem. These include perturbation theory (Fermi's golden rule), 
numerical integration of the time-dependent Schrodinger equation and semi- 
classical phase space analysis leading to stochastic ionization T, TF, l^l] ■ 
Rigorous approaches include Floquet theory and complex dilations [18, 19, |3T| . 
Despite this, there are few exact results available for the ionization of a bound 
particle by a realistic time-periodic electric field of dipole form E(t) ■ x (an 
AC-Stark field) for fields of arbitrary strength. The most realistic results we 
are aware of are based on complex scaling ([TSl [Tni [31]) and show ionization 
(for small electric field) of certain bound states of the Coulomb atom as well as 
defining resonances in some regions of the complex energy plane. 

The lack of rigorous results for large electric fields is true not only for re- 
alistic systems with Coulombic binding potential, but even for model systems 
with short range binding potentials [H HI |T7] . The most idealized version of the 
latter has an attractive (5-function potential in 1 dimension. The unperturbed 
Hamiltonian Hq = —d^ — 2S{x) has a bound state (t>o{x) = e~'^' with energy 
— 1, and explicitly known continuum states [8 . This model has been studied 
extensively in the literature, but the only rigorous results (known to us) con- 
cerning ionization involve short range external forcing potentials rather than 



dipole interaction; see however [31 [TH [53] for some rigorous bounds on the ion- 
ization probabiHty by a dipole potential for finite time pulses. Detailed results 
for compactly supported forcings were obtained in [71 [§[ [5S1 13] . In this paper 
we develop techniques to deal with physically realistic dipole interactions. 

We consider the time evolution of a particle in one dimension governed by 
the Schrodinger equation (in appropriate units): 

idt^Pix,t) = ( --^ -25{x)] ^p{x,t) + E{t)xi}{x,t) (1.1a) 



dx^ 

^ix,0) = Mx)eL^iR) (1.1b) 

Here, {x,t) e M x IR+, E{t) is real valued, smooth and periodic, E{t) = E{t + 
l-K Ilo). We prove that for E(t) a trigonometric polynomial, 

N 

E{t) = (Ene'"^* + E;:e-'"^*) , (1.2) 

n=l 

with N < oo the system always ionizes, i.e. for ip{x, t) solving ()1.1|) . 

lim / \ij{x,t)\'^dx = 0, VLeM+ (1.3) 

with the approach to zero at least as fast as 0{t^^) if ^o{x) E D L^(M). 
When E{t) is not a trigonometric polynomial (i.e. = oo in p.2p ). the Floquct 
Hamiltonian (see below) may have time-dependent bound states and ionization 
may fail. This is uncommon, but there are examples of time periodic operators 
where such bound states exist [3 125] . 

For general periodic E{t) and for compactly supported initial data, il){x^t) 
can be uniquely decomposed into a sum of (non-L-^) functions 

"Pi^'t) = J2 I]afcji'e-**''*$fc,,(x,t) + ^-,^2;,*) (1.4) 

fe=0 j=Q 

(with 1 < M < oo). The indexing in k is chosen so that Sct^ < '^ak+i- 
The functions $fcj(a;, t) are 27r/w periodic in time and exponentially growing in 
space, behaving near x — ±oo like the exponentially growing Green's function 
of the Floquet Hamiltonian. Note that Q&k < 0, otherwise the wavefunction 
would grow exponentially in time. If we define 

7fc = -Qo-fe, (1.5) 

then for small values of 7^, 2jk gives the dominant part of the ionization rate 
for the fc-th resonance. The smallest rate, 70, gives the overall ionization rate 
for most experimentally relevant times [5] ■ The ionization rate 70 will vanish as 
E(t) 0, while 7/c ^ 00 as E(t) ^ for A; > 1. 

The term "^Mixjt) is a remainder after collecting Af resonances (provided 
iplxjt) has M resonances), possessing an asymptotic expansion in time with 
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power law terms in for n > 1. ^'M(2:,i) is computed as the integral of 

a certain function around a branch cut; the resonant states ^k,j{x,t) are the 
residues of the poles of that same function (this M can not be greater than the 
number of poles). We note that the polynomially decaying component of the 
wavefunction, to which we shall refer as the dispersive part, has actually been 
observed experimentally [57j, although under significantly different^ physical 
conditions. 

Replacing E{t) by eE{t), (Tq and ^o,o{x,t) have convergent power series 
expansions in e when u;^^ ^ N. When e —> 0, we have e~'°"''*$o,o(a^, ^) — ^ 
e**e~l^l, the bound state of Hq and '$i{x,t) goes to the projection of il^{x,t) 
on the continuum states of i?o- This shows that the first resonance is the 
analytic continuation in e of the bound state. This rigorously justifies some 
standard physics calculations in [12l [16l [2l] (see also the forthcoming work 
[30j . from which we drew inspiration). The Fermi golden rule and multiphoton 
generalizations can be recovered in our formalism through perturbation theory. 

All other resonances must come from a = ~ioo. That is, as e ^ 0, 7^ — > 00, 
meaning that these states "ionize instantly" in the absence of the electric field. 
We conjecture that such states do not exist for this model. Indeed, in all other 
cases considered [315], such states do not exist, but our technique does not rule 
them out. See Remark l3.14l for more details on this point. 

1.1 Outline of the strategy 

Due to the fact that the binding potential d{x) has support {0}, the behavior of 
'ip{0^t) and the initial condition completely determine the behavior of the solu- 
tion. We shall therefore deal mostly with ip{0,t) which satisfies an autonomous 
equation, (|3.4p . That equation is well posed, and its solution can be extended 
to ip{x,t) solving (|l.ip . This is sketched in Appendix [Cl 

Our main tool is the study of the analytic structure of the Zak transform of 
ipiO^t) (with V(0,t) = for t < 0), 

Z[V'(0,.)](a,t) =^e'-(™/-)^(0,t + 2^jM (1.6) 

in the complex a domain. As will be made clear later, see Definition 13.11 and 
Proposition [3l2l it is sufficient to consider the strip < Sftu < lu, which we shall 
do henceforth. Unitary evolution of the wavefunction implies that Z[ip(0, ■)]{ct, t) 
is analytic in a for 3(T > 0. 

For E{t) = and (V'(x, 0)|e-l^l) 0, Z[^P{0, ■)]{a,t) has a pole at ak = 
— 1+ [l/uj\uj corresponding to the eigenvalue —1 oi Hq (with [l/wj the largest 
integer less than l/w). The residue at the pole is e~*Li/'^J'^*, jf -^^e consider a 
outside the strip < < uj, this pole is repeated at the points ak + muj (see 
Definition 13.11 in particular p.6cp ). 

^ In 1271 , the authors studied luminescent decay of dissolved organic materials after a pulsed 
laser excitation. 
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We show that when E{t) ^ 0, the poles give rise to the terms e~"^'=*$fcj (2:, t), 
with the residue at the poles corresponding (by a linear transformation) to the 
Fourier coefficients in time of ^k,jix,t). There is also a branch point at cr = 
which gives rise to the dispersive part of the remainder Mix,t). 

The proof of complete ionization, (|1.3p . involves proving that there do not 
exist any Floquet bound state (nonzero solution to (jl.Sp V This is done by 
solving the Schrodinger equation without the (5-function at zero, and showing 
that solutions which decay as x — > —00 can not be matched continuously at 
X = to solutions which decay exponentially as x +00, implying that Qcr < 0. 

1.2 Statement of results 

We consider the Schrodinger equation with a time periodic Stark Hamiltonian 
(jl.ip on ]R^+^. E{t) is continuous, and is given by (|1.2p . We prove two theorems: 

Theorem 1 (Ionization) Suppose E{t) is a trigonometric polynomial, i.e. 
En = for n > N (with N < 00). Then for any ipo{x) € L^(M) ionization 
occurs in the sense of U.!^) . If ipai^) G i^(R) n i^(M), then the approach to 
zero is at least as fast as t^^ . 

A key tool in proving Theorem [T] is the following result on the structure of 
ipixjt). This result holds even if = cx) in (|1.2p . 

Theorem 2 Suppose 'ipoix) is compactly supported and in (finite kinetic 
energy). Then, the solution ip{x,t) of (|l.ip . can be decomposed as: 

M rifc 

V'(x,t) =^^afc,jt^e-^*'=*$fc,,(x,i) + *M(a:,t) (1.7) 

fc=0 j=0 

where 5(Tfc < 0, $fe.j(x,i) is lixjio periodic in time and continuous in x. Here 
M is finite, but the number of resonances may be infinite. The resonant term 
^k,j{x,t) is a Gamow vector and grows like the exponentially large Green's func- 
tion as \x\ — > 00. ^k.j{x,t) is continuous in x and t, differentiable in t and x 
except at X — Q. For j = 0, $fc.o(x, t) is an eigenvector of the Floquet Hamilto- 
nian: 

(^~idt ~ ^ - 2^(x) + E{t)x^ $fe,o(x, t) = ak^kfiix, t) (1.8a) 

^kA^,t) = ^k,o{x,t + 2TT/iu) (1.8b) 
The remainder \l/jv/(a;,i) has the following asymptotic expansion in time: 

00 

^m{x, t)^Yl E Dj,n{x)t-"/^ (1.9) 
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Finally, when G (0,w) (i.e. — Qo'fc —0), then — and ^kfli^it) is 
an L'^ -eigenvector of the Floquet Hamiltonian and decays with x, 

lim $fc,o(a;,0 = lim $fe,o(a;,t) = (1.10) 

X — — OO X — ^C30 

In this case, the functions ^k,Q{x,t) and M{x,t) are orthogonal. 

The Gamow vectors described in Theorem [2] are time dependent, and can 
be written expUcitly as follows: 



f I L i^&k+nujx —inLjt —ib(t)x—ia(t) <' fl 

$j.,fxt) = <^ , - fill) 

"''J ^ ' ' 1^ ^ ^R^—i^dk+nuix ^ — inujt ^-~ib(t)x~ia(t) 2; > Q 

The coefRcients V'n'^ vary with j and k. The functions a{t) and &(i) are defined 
in p.l2p . It should be noted that if 3(7^ = 0, then V'n = for n > (otherwise 
^k,Q{x,t) would not be in L^, and (|1.10p would be violated). 

Remark 1.1 The PDE (|1.8a[) is formally overdetermined if = 0, since it 
has three boundary conditions ( (jl.lO|) and (jl.Sbp ). This makes nonzero solutions 
to (jl.Sp unlikely, although there may be some special forms of E{t) for which 
such a solution can be found. The proof of Theorem [1] is essentially a proof that 
in the case of E{t) a trigonometric polynomial, there are no nonzero solutions 
in this setting. This also implies that 7^ > for all k 

Remark 1.2 Although Theorem [5] applies only to compactly supported initial 
conditions in , if all compactly supported initial conditions ionize, then all 
initial conditions in L^(]R) will ionize. This is a simple application of the fol- 
lowing well known resultd to the operator family T{t) — \^_]^ L]{x)U{t) {Ult) is 
the propagator for ()1.1|) ): 

If T{t) is a uniformly bounded family of bounded operators on i^(R), and if 
T{t)u ^ for M in a dense subset of i^(M), then T{t)u for aU u € L'^(R). 

Remark 1.3 11-00(2^) is not compactly supported, then a similar decomposition 
can be provided, but with extra terms. Essentially, if Z[e'^^'=*ipn{x)]{<T,x,t) is a 
ramified analytic function of a, each singularity will lead to a similar singularity 
in Z[ip{x,t)]{a, x,t), leading to other exponential decay terms. Resonances are 
poles in a, the location of which does not depend on the initial condition. 

Remark 1.4 We believe the dispersive part 'i>M{x,t) is Borel summable, al- 
though this does not follow from our results. To show this, one needs to find 
exponential bounds on Z['ip{0,t)]{a,t) as Gcr -co, which would also show 
that there is only one pole, ctq, the analytic continuation of the bound state. 



A simple proof: for any ip in L'^{R), find a 11 a distance e away from i/"- Then ||r(t)?/'|| < 
\\T{t){ij - u)\\ + \\T{t)u\\. Note that ||r(t)(u - < \\T{t)\\ e < Ce (with C the uniform 
bound on ||r(t)||), and for sufficiently large t, ||T(t)M|| < e as well. Thus ||T(t)V'|| < (C + l)e 
for sufficiently large t. 
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1.3 Equivalent formulations 

Here we describe some equivalent formulations of This material is essen- 

tially taken from chapter 7 of [TU]. We will use (|1.14p in the proof of Theorem [1] 
and (|1.13p in the proof of Theorem We first define some auxiliary functions: 

a{t) = [ b{sfds = aot + a^{t) (1.12a) 



b{t) = y ( J^e™"* + _^e'™"* ) (1.12b) 
^-^ \ into —inuj ' 



n=l 



°° / F F 
c(t) = 2 y — ^e™"* + , . " ,, e-'""* = y (C„e""* + C„e-™'^*) 



(1.12c) 

where av{t) is 2tt/uj periodic and has mean 0, and ao = (u;/27r) bis)^ds. 
Note that (l/2)c"(i) = b'{t) = E{t). 

Define il>v{x,t) = e+'-^it) e+^-^m^-cit))^^^ „ c(<),i); then the following equa- 
tion for ijjv is equivalent to 

idtiJvix, t) ^ (^--^ - 2(5(x - cit))^ 4>4x, t) (1.13) 

This is the velocity gauge, and the equivalence can be verified by a computatiord. 
Similarly, there is an equivalent equation in the magnetic gauge. We obtain it 
by setting i}^B{x,t) = e+'''(*)e+*''(*)="V(a;, <): 

idt^B{x,t)= {-^-25{x) + 2ib{t)d^^B{x,t) (1.14) 

Remark 1.5 Suppose that either ipB{x,t) or ipy{x,t) are time-periodic solu- 
tions of (|1.14[) or (|1.13p . Then tp{x,t) is a time quasi-periodic solution of (|l.ip . 
and e'^°'"*ip{x, t) is time-periodic. 

Of course, all this follows only after showing that (|1.1|) . (I1.13P or (|1.14p are 
well posed. This is discussed in Appendix [C| the basic idea is to solve Duhamel's 
equation for ?/'(0, t), namely (|3.4p . and then extend the solution to all x. 



1.4 Organization of the paper 

In Section[2l we assume Theorem[2]to be true and use it to prove Theorem[T] In 
Section [3] we prove Theorem [21 In Section 3] we make some concluding remarks, 
and discuss possible directions of future research. Some technical material is 
presented in the appendices. 

■^Equation HI. 131 1 differs from wliat one finds in 1101 . In |10| . the authors take 6(t) = 
Jq E{s)ds and c(t) = Jj* b{t)dt, which imply that c(t) = c{t) + cq + Cvt. This does not change 
the essential feature that {l/2)c"(t) = b'{t) = E{t). 
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2 Ionization 



Based on Theorem [2l we will to show that the Floquet equation (jl.Sp in the 
magnetic gauge has no nonzero solutions with Qcr — which satisfy (jl.lOp . This 
implies ionization for compactly supported initial data, which by Remark 11.21 
implies ionization for all ipo{x) € 

In Section [^m we solve (|1.14p without a binding potential (the —25(x) term) 
and characterize the solutions. We then assume a bound state ^k,o{x,t) exists, 
expand it in an appropriate basis, and derive necessary conditions on the coef- 
ficients to meet the boundary conditions (decay at x = ±oo and continuity at 
X = 0). 

In Section 12.21 we use the characterization of solutions we constructed in 
Section [2. II and show for E(t) a trigonometric polynomial that there are no con- 
tinuous, nonzero solutions to (|1.14p which vanish at x = ±oo. The basic tech- 
nique is to analytically continue, in the t variable, both tpB{0-,t) and iPb{Q+, t) 
(which must coincide) and use the Phragmen-Lindelof theorem to show that an 
associated function must be entire and bounded (and therefore constant). This 
implies that any localized solution to p.Sp is zero, and ionization occurs. 

2.1 Solutions to the free problem 

By Theorem [2l we need to show that (|1.8p has no nontrivial solutions. In the 
magnetic gauge, this is the same as showing that if ^k,o{x,t) solves 

<7k^kAx,t) = i-idt ~ dl ^ 25{x) + 2tb{t)d^)^k,o{x,t) (2.1) 

with boundary conditions (ll.lOp and (jl.Sbp . then $^^0(2^,^) = 0. 

We begin by solving (|2.ip without the (5-function binding potential (and 
letting cr = (Tfc, which causes no confusion in this section), 

a^P{x, t) = {-idt -dl + 2ib{t)d^)^j{x, t) (2.2) 
Taking i/j{x,t) = e^^(f\{t) as an ansatz, we obtain an ODE for f\{t): 

dtvx{t) = ~i {-a -\^ + 2iXb{t)) ^x{t) (2.3) 
This has the following family of solutions (recalling that c'{t) — 2b{t)): 

ip^{t) = e-'^^'e^<'^ (2.4) 

To ensure 2tt/uj periodicity in time, we must have (— <t— A^) = mco, m G Z. This 
implies that A = zbi-^/mw -|- a (with the branch cut of ^/z taken to be — iM+). 
Therefore, (|2.2p has the family of solutions: 

fmM^^t) = e±A™:.g-«m^tg±A„c(t) (2.5a) 

\„i = —i\/(T + muj (2.5b) 
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2.2 Matching solutions 

Given the family of solutions to (|2.2p . we can attempt to solve (|1.14p . Applying 
Theorem [5J we have three boundary conditions to satisfy: 

$fc,o(0,i) = $m(0-><) = 3'fe,o(0+,i) (2.6a) 

5^$fc,o(0+,t)-6),$fc,o(0_,<) = -2$fe,o(0,i) (2.6b) 

lim $fc,o(-2;,t) = lim $fe,o(+x,i)=0 (2.6c) 



a: — ^oo X — ^oo 



Consider now a solution ^k,o{x,t). We can expand (formally) iMx,t) in 
terms of the functions (pm,± in the regions a; < and x > separateljO: 

a, .^_ f 12mezi^^,+'Pm,+ {x,t) ip^,^{x,t)), x<0 



Emez(V'm, + ¥'m,+ (a;,i) +^m,_</?m,-(2;,i)), X > 

For TO > (recalling ak G [0,'^.') and examining (|2.5bp '). the functions 
ipm,±{x,t) are oscillatory in x as a; — > ±00. Thus, if the coefficients ipm ± 
(jn > 0) were not zero, then $^^0(2;,^) would not decay as a; — > ±cx), violating 
(IM . 

Similarly, we observe that ipm,+ {x, t) are exponentially growing when to < 
as X — !■ +00, so V'm + must similarly be zero. The same argument applied to 
the region x < shows that ■(/'m + must be zero when to < 0. Therefore after 
dropping the ± in the coefficients V'm ±j we obtain the result we seek. 

Thus, we find that we can actually write $fe_o(a^, t) as: 

$,.o(-,^) = ( &"<"^l''™^i"'S' "tl (2.8) 
I Em<oV',„¥'m,-(a;,t), x>0 

with both sequences "0™^ m P. Although this derivation is purely formal, it is 
proved in Appendix IB] It also motivates (|l.lip . 

Substituting ()2.8p into the continuity condition p.6ap yields: 

V'™e"""'^*e^'""(*) = V'me-'""'^*e-^'"'=(*) (2.9) 

m<0 m<0 

Proposition 2.1 Suppose E{t) is a trigonometric polynomial with highest mode 
N, that is E{t) = E!Li(^«e'""* + i?«e-'"'^*). Set z = e"*'^*. Then $fc,o(0,i) 
has the decomposition: 

^kA^.t)^ f{z)+g{z-^) (2.10) 

The Junctions /(.) and <?(•) are entire Junctions oj exponential order 2N , and 
(7(0) = 0. This shows in particular that $/c,o(0,t) is continuous. 

The correspondence between $fc^o(0, i), f{z) and g{z) is as Jollows. Let 
denote the j'th Fourier coefjicient oJ ^k,o{0,t), that is $fc,o(0,t) = V'je*.'"*. 
Then letting fj, gj be the Taylor coejficients oJ J{z), g{z), we Jind fj — ip-j Jor 
j > and gj = ipj Jor j < 0- 

■^The validity of the expansion is proved in Lemma IB. 2l in Appendix [B] 
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The proof of this fact uses results from Section [31 and is deferred to Appendix 
[XI Finally, we state a result we use, proved in most complex analysis textbooks, 
e.g. [IS]. 

Theorem 3 (Phragmen-Lindelof ) Let f{z) be an analytic function of ex- 
ponential order 2N , that is 1/(2) | < Ce'^ . Let S be a sector of opening 
smaller than tt/2N . Then: 

sup 1/(^)1 > sup |/(z)| 

z€dS zGS 

We are now prepared to prove the main result. 
Proof of Theorem [7[ 

We describe first the case iV = 1 now (i.e. E{t) — Ecos{ujt); the case of 
arbitrary N is treated below) . The key idea is that we can use (12. 8p to obtain an 
asymptotic expansion of $fe,o(0+, t) and <&/£.o(0_, <) in the open right and left half 
planes in the variable z = g"*'^* (respectively); to leading order ^i^ Q{0_,t) ^ 
^/'^z'"e-^l^^l and $fe,o(0+,t) ~ Vm^'^e"^!^^! (note that to and C may be 
different). This asymptotic expansion shows that f{z) decays exponentially 
along any ray z — re"^ in the open left or right half planes. 

In fact, the asymptotic expansion allows us to observe that f{z) (the part 
of ^k.oiO, t) which is analytic in z) must be boimded except possibly on the line 
iM. Theorem [3] combined with Proposition 12.11 allow us to conclude that f{z) is 
bounded on the line zM. This shows f{z) is bounded on C and hence zero. 

Since f{z) is zero, $fe,o(0,t) = g{z) ^ Omz^^'' for some Af e N (since g{z) 
is analytic). But we previously showed also that $fe o(0,t) ~ V'm-^™^"*^'^^' • 
Two asymptotic expansions must agree to leading order; the only way this can 
happen is if g{z) = $fc,o(0,<) = 0. 

The main difference between the case iV = 1 (monochromatic field) and N > 
1 (polychromatic field) is that instead of the exponential asymptotic expansions 
being valid in the left and right half planes, they are valid in sectors of opening 
tt/N; to show this we need to apply Theorem [3] to the boundaries of these 
sectors. 

We now go through the details. 
Step 1: Setup 

Let ^k,o{x,t) be a solution to (|1.8p . By the hypothesis of Theorem [l] we 
let E{t) be a nonzero trigonometric polynomial of order N. Let z = e~*'^*. Let 
e^(^) = {Cjz3 + CjZ-^) where the Q are the coefficients from (|1.12cp . 

We apply Proposition l2.1l to <i>fc o(0, i) and (|2.9p to obtain: 

<^>kAO,t) = f{z)+g{z~') 

= Y: ^^.z"e+^-^(^' = (2.11) 

r?i<0 T?i<0 

The first equality holds by (|2.10p . the second by (|2.8p with a; = 0. A priori, 
equality holds only when |z| = 1. However, both of the latter two sums are 
analytic in any neighborhood of the unit circle in which they are uniformly 
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convergent. Thus, f{z) + g{z~^) is the analytic continuation of the sum if the 
sum is convergent in some neighborhood containing part of the unit disk. 

For the rest of this proof, we make the following convention. The functions 

ip'"'^{z) are defined by 

V'^(z) = Yl ^™^™e+^-«(^) (2.12a) 

m<0 

V-^Cz) = Yl V-^^^e-^-^^^) (2.12b) 

m<0 

for those z for which the sum is convergent. 
Step 2: Convergence of the sum 

We show now that the sum in (j2.11|l is convergent in a sufficiently large 
region. 

For \z\ > 1 and ^€{z) > 0, consider the sum X)r?x<o ^m-^'"^"^'"'^*^^''- 
this region, since 3fJ£(z) > 0, we find that e^"'^^^^ < 1. The coefficients V'm^ 
are bounded uniformly in m (since they form an P sequence). For \z\ > 1, 
z™ is geometrically decaying as m —oo. Therefore the series is absolutely 
convergent when \z\ > 1 and 5R£(z) > 0. 

The same statement holds with '^^^Q^mZ"^e^'^'^'^^^'> in the region where 
m{z) < 0. 

Let us define the following sets: 

5*+ = Connected component of 5Mn {z e C : \z\ > l,^€{z) > 0} 

Connected component of 5Mn {z e C : |z| > l,^(L{z) < 0} 

A plot indicating the structure of these sectors (for a particular choice of £(z)) 
is shown in Figure 1 for the case where N = 2. 

By Proposition l2.11 we see that ip-^{z) is analytic in 5*+ and ip^{z) is analytic 
in S~ , since the sum in (|2.12p is convergent there. 

We now show that 5*+ and must be unbounded since £(z) is not constant. 
First, note that £(z) — £(z~i). As in the Schwarz reflection principle, define 
B = S+U (5+)-i. Clearly, m{z) = for z G Ob- If S+ is bounded, then B is 
bounded as well. By the real max modulus principle, 5R£(z) must be zero inside 
B, and hence 5R£(z) is bounded everywhere, which is impossible. 

Finally we show that the regions S'^ and S~ "fill out" to open sectors as 
|z| — > oo. That is to say, if S is some sector in which 5Rz^ > 0, then for any ray 
{re'^ : r > 1} contained in S, there exists R — R{9) so that the truncated ray 
{re'" : r > R{9)} C S+ . 

Without loss of generalitj0, let us suppose that Cn G M+. For very large 
|z|, we write €{z) = J^^^i CjZ^ + Cjz-^ = C^vz^ + O(z^-i). Then setting 
z = re''', we find that r'^^ire''') = Cnc''^'' + 0{r-^). Thus, for r sufficiently 
large and N9 ^ (2m+l)7r/2, we find that r~^£(re'^) has either strictly positive 

^ Suppose Cjv = pe'^ . Then rather than choosing z = e"^^, we would substitute z = 
^i(^t-e/N)^ 
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real part or strictly negative real part. In particular, if \N9^'k/2\ > e, then 
there exists an R = R{e, 6) so that 3?r~^£(re'^) is bounded strictly away from 
zero. 

Motivated by the above, we define the following subsets of C (with j = 
0...7V-1): 

A+^ = {re'^ ■.r>R{t,e), 

9 e [-7r/2N + 27rj/N + e, n/2N + 2nj/N - e]} (2.13a) 



A-^ = {re'' ■.r>R{e,e), 

9 e [-n/2N + 2n{j + l/2)/N + e,7r/2iV + 27r(j + l/2)/A^ - e]} (2.13b) 

Clearly, for sufficiently large R, A^^ \ Br C 5+ and A'J ^ \ Br C . Here, Br 
is the ball of radius R about z = 0. 
Step 3: Asymptotics of f{z) 

We now show that f{z) — 0. We begin by writing f{z) as follows: 





oo 


oo 

- E 5"-"" - 


m<0 


/(^) 


oo 
n=0 


oo 

~E5n^""- 
n=l 


^Ev^. 

m<0 



,^z"e-^"'^(^),z G 5+ (2.14a) 



i^mg+A„e(.)^^^ 5.- (2.14b) 

We let 5fc, fc = 0, . . . , 27V + 1 be a set of sectors of opening tt/{2N + 1) 
arranged in such a way that the boundaries of Sk avoid the rays re*^^^^+"'^'/^^. 
Therefore, for sufficiently large \z\, the boundaries of Sk are contained in either 
or A~^ except for a compact region. On dSk, f{z) is decaying as \z\ 00, 
by a simple examination of (|2.14p . Since f{z) is entire (unlike ip{z)), f{z) is 
also bounded on dSk even for small z. 

We have shown that f{z) is bounded on dSk- Applying the Phragmen- 
Lindelof theorem, f{z) is therefore bounded on Sk- Since U^^^Sk = C, we find 
f{z) is constant. Since we know that along any ray contained in ^^g, fiz) is 
decreasing, we know /(z) = 0. 

Step 4- Asymptotics of g{z) 

We now show that g{z) = 0. We rewrite (|2.14[l with g{z) on the left side. 

00 

Effn^-" = E V'^^'"e-^'"«(^),ze 5+ (2.15a) 

n— 1 m<0 

00 

J29nz-" = E V'™^'"e+^'"^(^),z e S- (2.15b) 

n— 1 m<0 

Since the left sides of (|2.15aP and (j2.15bp are (convergent) asymptotic power 
series (for sufficiently large |z|), while the right sides of (|2.15ap and (I2.15bp are 
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(convergent) asymptotic series of exponentials, we find that the right side decays 
much faster than the left side. This is impossible unless both sides are zero. 

□ 



3 The Floquet Formulation 

In this section we prove Theorem [2] To do so we define an auxiliary function 
Y{t) = %p{c{t), t) and derive a closed integral equation of Volterra type for it via 
Duhamel's formula. We then apply the Zak transform in time to the integral 
equation for Y{t). This yields an integral equation of compact Fredholm type 
for Z[y]((7, t), the Zak transform of y{t). The integral operator is shown to 
be analytic in a. Applying the analytic Fredholm alternative to this equation, 
shows that Z\Y][a,t) is meromorphic in a^/"^. The poles corresponds to reso- 
nances or bound states, while the branch point corresponds to the dispersive 
part of the solution. 

In Section 13.31 we extend these results from a; = to the entire real line. 
We show that the wavefunction, considered in the magnetic gauge, can be de- 
composed in the form (|1.7p . If So-fc = 0, then IRdk € (0,a;) and <^k,o{x,t) 
corresponds to a Floquet bound state. The remainder "^M{x^t) decays with 
time, in particular ||\l/M(a;, <)|j2,cx, < 0(t~^^^) as t ^ oo. 

3.1 Setting up the problem 

Here we work in the velocity gauge and study (|1.13p . Recall that c{t) is 27r/w 
periodic. We rewrite (jl.ip in Duhamel form, using the standard Green's function 
for the free Schrodinger equation, (47rzt)~^/^e" 

+ 2i f [ e^v{ '^^,7'''fA s{x'^c{t'))Mx',t')dx' '^f (3.1) 
Jo Jr \ 4(t-r) J y/ATTi{t-t') 

where we have defined: 

i^,^oix,t) ^ e'^'*Mx,0) = f (47rit)-i/2e*h-"T/^V«(^',0)da:' 

Jr 

Computing the x' integral explicitly and changing variables to s = t — t', we 
find: 

+ 2i/exp ^ \ ^U^c t-s ,t-sH^ 3.2 
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We now substitute x — c(t), to obtain a closed equation for ^j^{c{t), t): 

IT /■* f i(c(t) - c(t- s))^\ , , , , ^ds , , 
Setting yo(i) = '0t),o(c(i), and = ip{c{t),t) we obtain: 

r(0 = VMO + exp r« - ,3.4, 

Here the derivation of (|3.4p is formal; a sketch of a rigorous derivation can 
be found in Appendix [C] The basic idea is to solve (j3.4p and then use ()3.2p 
to extend the solution to all a;. Showing the solution is in for each t is 
accomplished by stationary phase, see Appendix [Cl 

The main tool of our analysis will be the Zak transform. 

Definition 3.1 Let f{t) = for t < and \f{t)\ < Ce"* (a G R+). Then 
f{t) is said to be Zak transformable. The Zak transform of f(t) is defined (for 
SJfj > a) by: 

Z[f]M = ^e''^(*+2-^»/(t + 27rj7c^) (3.5) 

and by the analytic continuation of (|3.5p when 3(t < a, provided that the ana- 
lytic continuation exists (treating Z[f\{a,t) as a function of a taking values in 
L2([0,2Vc.],dt);. 

Proposition 3.2 Z[f]{a,t) has the following properties: 

fit) = Lo-^ e-''''Z[f]{a,t)da (3.6a) 

J 1/3 

If Z[f]{a,t) is singular for 3cr = /?, this integral is interpreted as the limit of 
integrals over the contours [i{(3 + e),i{l3 + e) + uj] as e — )■ from above. 

Z[f]{a,t + 27r/Lo) = Z[f]{a,t) (3.6b) 

Z[f]{a + LO,t)^e'^'Z[f]{a,t) (3.6c) 
If p{t) is 2'k/uj -periodic, then: 

Z[pf]{<j,t) = p{t)Z[f]{a,t) (3.6d) 

With the exception of (|3.6ap , these results all follow immediately from (|3.5p . 
See Remark 13.51 for an explanation of (|3.6ap . 

Remark 3.3 Suppose f{t) is Zak transformable, and uniformly bounded in 
time {a = 0). Suppose further that the analytic continuation of Z[f]{a,t) has 
a singularity (say at cr = 0). Then (|3.6cP still holds, in the sense that for any 
direction 0, Z[f]{cr + lo + Oe'" ,t) = e"^*Z[f]{a + Oe''^ ,t). 
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Remark 3.4 More information on the Zak transform can be found in, e.g., [TTJ 
p.p. 109-110]. Our definition differs slightly from that in [Til by allowing a to 
take complex values. 

Remark 3.5 One can relate the Zak and Fourier transforms as follows. Let 
f{k) = j e^'^* f{t)dt denote the Fourier transform of f{t). Then: 

= |i^/(a + na;)e— * (3.7) 

riGZ 

The Poisson summation formula, applied to (j3.5p . yields (|3.7[) . Eq. (|3.6aP fol- 
lows immediately from p.7p . This relation implies that our approach is equiv- 
alent to the Fourier/Laplace transform analysis done in [6l [71 [15]. The Zak 
transform is used simply for algebraic convenience. 

We proceed as follows. Applying the Zak transform to (13. 4p yields an integral 
equation of the form 

y[a,t)^yo{a,t)+K{a)v{<j,t) (3.8) 

with y{(T,t) — Z[Y]{a,t), yQ(a,t) = Z[Yo]{a,t) and K{a) the Zak transform of 
the integral operator in (|3.4p . K{a) will be shown to be meromorphic in <j as 
a compact operator family from L'^{S^,dt) — > L'^{S^,dt), except for a branch 
point at cr = 0. 

We then use the Fredholm alternative theorem to invert (1 — K{a)). Once 
this is done, we find: 

y{a,t) = il-K{a))-'yoi<J,t) (3.9) 

The poles of {1 — K{a))^^ correspond to resonances, and a branch point at cr = 
corresponds to the dispersive part of the solution, i.e. the part with polynomial 
decay in t~^/^. 

To begin, we determine the analyticity properties of Z[yo](c, t). 

Proposition 3.6 Suppose ipoix) is smooth and compactly supported. Then near 
a — 0, yo{<J,t) has the expansion: 

Z[Yo]ia,t) ^ yo{a,t) = ^a'^/^ f Mx)dx + f{<J^'\t) (3.10) 

The function f(a^^^,t) is analytic in cr^^^, and is in L^{S^,dt). Also, for some 
constants Ci and C2, we have 

Here, is the set [0, 27r/w] with periodic boundaries (so that = 2t\ juj in S]^). 
In fact, the same conclusion follows for Z[il>(x + c{t), t)]{a, t). 
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Proof. Consider Yo{t) = Vi),o(c(i), i) for t > only (and Yo{t) = for i < 0). 
Then (with slight abuse of notation) we write 



JR 

Computing the Zak transform yields: 



jez 
s^'=^(*)i/.o(fc) 



dk 



ikc{t) 



E 

Ln6Z 



(1/2) / My)dy + <y^'^'(i/2) / (e-^l^(*)-^i - i)V^o(y)dy 



E 



2\/c + 



^V^^+^W(t)~y\^^(^y)dy (3.11) 



The interchange of the sum and integral between lines 1 and 2 is justified (for 
Scr > and t fixed) since the sum over j is absolutely convergent, as is the 
integral over k. The result is valid for arbitrary a by analytic continuation. 

The change inside the square brackets between lines 2 and 3 comes from 
the Poisson summation formula in the t variable, and the fact that the Fourier 
transform of xr+ it)e'^^^+'''>* is -i{k'^ + a + C)"^(with C dual to t). 

The first term on the right side of ()3.1ip agrees with that in (I3.10|) . Since 
(e~^l'^*^*''~^l — 1) is analytic in a^^'^, the second term is analytic in u^/^. The 
second and third (which is analytic in a) terms become f{a,t). 

Since ipo{x) is supported on a compact region, \c{t) — y\ is bounded (say by 
C2) and exponential growth follows. 

This result follows for all x rather than simply simply by translation in- 
variance of 6*^==*. □ 



We now determine the Zak transform of the integral operator in (|3.4|) and 
compute the resolvent of it. 

3.2 Construction of the resolvent 

We now apply the Zak transform to p.4p to construct an equivalent integral 
equation. 

Proposition 3.7 Let f{t) be Zak transformable. Consider the integral operator: 



Kvf{t) 



exp I 



{c{t ) - cjt - s)f 
is 



fit~s) 



ds 

7^ 



(3.12) 
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Then if > 0, we find: 
Z[Kvf]iaJ)=K{a)f{a,t) 



TT Jq \ As J 

Proof. We rewrite p.l2p as: 

Applying Z to both sides of (13.14^ yields 
Z[Kvf] (rx, t)^Yl ^'^^'^^^'1^^ {Kvf\ it) 

= {IY.--^^-'^''^^ I exp (.^^^^^^^^) fit + 2nj/. - s)x.. (.)^ 



ds 

XWL+(.S)—J= 



, expf.M^^(i^)e-^Z[/](.,i-.)i^ (3.15) 

This is what we wanted to show. □ 

We now show that the operator K{a), constructed above, is compact. We 
decompose K{a) as Kp{a) + K^icr) (defined shortly), and treat each piece 
separately. 

Proposition 3.8 Define Kf{(j) : L'^{S^,dt) L'^{S^,dt) by: 

Then, Kp{a) is compact and analytic for Scr > 0. It can be analytically con- 
tinued to 3(T < 0, (T ^ 0, and the continuation has a a~^l'^ branch point at 
cr = 0. 

Proof. We compute this exactly by expanding f{t) in Fourier series and 
interchanging the order of summation and integration: 
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This is valid for Qcr > 0, as well as 5cr — but in this case we must treat the 
integral as improper. 

Thus, in the basis e~'""*, this operator is diagonal multiplication by {a + 
nw)"^/^. Compactness follows since the diagonal elements decay in both direc- 
tions. Analyticity for cr ^ follows by inspection of the right side of (|3.16p . and 
choosing the branch cut of \J<j + nio to lie on the negative real line. □ 



Proposition 3.9 Define Kl[(t) : L^{S^,dt) L^{S^,dt) as: 



KL{a)f{t) 



exp I 



{c{t) - c{t - s)f 



4s 



ds 



(3.17) 



Then Kl{<j) is compact for 3(t > and analytic for 3(t > 0. It has continuous 
limiting values at 3(t = 0. 



Proof. We rewrite p.l7p as: 



2Tr/u 



E 

fc=0 



exp I 



{ c{t)-c{t-s)y 
4(s + 27rfc/w) 



^ia{s-\-27rk/Lj) 

\/ S + 2Trk/uj 



f{t-s)ds (3.18) 



Provided that 5cr > 0, the sum is decaying at least as fast as k~^^^. Each term 
in the sum is continuous. Thus the sum is absolutely convergent to a smooth 
function in t and s, which is analytic in a (thus the limit is analytic except 
possibly when 3(t = 0). The region of integration is compact, and so is K^^a). 

□ 

We now analytically continue Kl{<t) to the strip < "Sia < lo. 
Proposition 3.10 Let K' {a) he the integral operator defined by: 



K'i<T)f{t) 



1 



exp 



k'^{t,s)f{t-s)ds 

{c{t)-c{t-s)f 
4p 



1 



(3.19a) 



^ (3.19b) 



dp 



where C is a contour along the real line in the upper half plane which avoids the 
singularities of the integrand at p = and p = i{s + 2T:n/uj) (see the proof for 
a specific example). 

Then K'(a) is analytic valued operator on < SRcr < lo, and vanishes as 
$5(7 — !■ oo. For every a, K'{a) is compact. Furthermore, K'{a) is the analytic 
continuation of (cr) . 

Finally, for a — ~iX (with SRA < 0) or a — —iX + uj, K(a) is analytic in the 
parameter X^/"^ or (a — uj)^/"^ . 

Proof. 

Step 1: Analyticity 
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When writing the contour of integration as K + Oi, we actually mean the 
integral along some contour in the upper half plane which avoids the singularities 
(at p = and p — i[s + {27r/uj)n]) of the integrand but remains close to R. 

For instance, let 7fl(i) = tujR/2n for t g R \ [— 27r/w, 27r/w], and 7r(<) = 
Re'^'^-^'^^+'^'^^/'^'i for t e [-27r/a;, 27r/w]. That is, 7fl(t) travels along the real line, 
and circles upward around the disk of radius R. The integral is then defined as 
limi?^o /^^ ■ dp. 

To compute the behavior of the integral, simply take R ~ I-k jio: 



B+Oi 



1 - e'^P' 



exp 



{c{t) ~ c{t ~ s)f 



1 - e'^P' 



exp 



Ap 

{c{t)-c{t-s)f 



dp 



- 1 



2m 



exp 



4p 

{c{t)-c{t-s)Y 
Ais 



dp 



- 1 



(3.20) 

VIS 



The integrand in the first term is analytic since p stays away from (thus 
avoiding the essential singularity at p = 0). It is exponentially decaying both 
for large positive p (at the rate e*^"'^'^-'^) and for large negative p (at the rate 
g-crpy If gfj^ = or = uj, the integrand still decays at the rate p^^l"^ , which 
is integrable. 

The last term is singular, but integrable at s = 0, and analytic elsewhere. 
Thus, k'^{t,s) has only a singularity of order s~^/^, and is analytic elsewhere. 
This shows that K'{a) is a compact family of operators, analytic on a. 

Step 2: Vanishing of the operator as 3(t — > +oo 

We examine p.20p . The first term vanishes as 5(t — > oo by the Riemann- 
Lebesgue lemma. The second term vanishes since e^'^'' does. Thus, k'^{t, s) — > 0, 
and so does K'{a). 

Step 3: Continuation of K{a) 

To show that K'{a) = K{a) if Scr > 0, we simply move the contour of 
integration in (j3.19b|) upward and collect residues: 



1 - e"P- 



, {c{t)-c{t-s)r 
exp I ; I - 1 



4p 



lim 



B+i27rAr/i. 



1 - e"P- 



exp 



dp 

VP 

{c{t)-c{t-s)r 



Ap 



dp 

VP 



N 



2iTi 



T(6*+27rj/u;) 



J=0 



exp 



(e(t)-e(t-s))- 
Ai{s + 2-Kjlui) 



1 



2-Ki 



,m(s+27rj/aj) 



J=0 



cxpi(-(^j--(^7)f 1-1 
^ ' Ai{s + 2-Kjlw) 



^i{s + 27rj7w) 
1 



i/i(s + 27rj7a; 
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We then integrate this kernel against an L'^{S^, dt) function f{t) and obtain: 



2n/u 



I LO 
TT 2111 

exp 



2-Ki 



^ia{s-\-27rj /uj) 



J=0 

Ai{s + 2-Kjluj) 
exp 



- 1 



^i{s + 2-Kjluj) 
4s 



f{t~s)ds 



e'^'fit - s) 



ds 



This is in agreement with (|3.17p . Hence, K'{a) = Kl{(j) for Scr > 0, SRcr G (0, uj) 
and therefore K'{a) is the analytic continuation of Kl{(j). 
Step 4- Singularity at a — 0,uj 

We now wish to show that K'{'-iX) is analytic in ^/X for a = —iX, and 
similarly that K{—iX + u!) is analytic in A. To do this, we proceed as in Step 3, 
but push the contour down instead of up. We rotate the contour 71 U 72 U 73, 
with 71 = [—ioo — R, — i?], 72 which goes around the unit circle of radius R in 
the upper half plane (as in step 1), and 73 which is [R,R — ioo]. This lets us 
avoid concerning ourselves with the singularities of the integrand; the important 
behavior is the decay near p = —ioo. 

Note that the integral kernel of K'{-~iX) is given by 



1 - e"P- 



/71U72U73 

while that oi K' {—iX + uj) is given by: 



exp 



k' 



{t,s) 



71U72U73 



^LUp—lLj; 



exp 



{c{t) - c{t - s)y 

Ap 



m-cit-s)f 

4p 



dp 

VP 



dp 



First, observe that the integral over 72 is analytic in A, provided R ^ uss. Thus, 
choosing a different R for cos < (3/4)7r and ujs > (l/4)7r shows analyticity in A. 

We consider the case a — —iX, the case a = —iX + uj being treated similarly. 
We now observe that, for ^p — R (the same argument applies to dip — —R), 
the integrand (over 73 or 71) becomes a Laplace transform: 



-iXp 



73 



— ioo 



exp 



(c(t)-c(t-s)) 



pUj{p+Fl) — iujs 



4p 
exp 



1 



dp 

VP 

{c{t)-c{t-s)r 



4{p + R) 



dp 



VpTR 



(3.21) 



We then observe that we can rewrite 



exp 



' icjt) - c{t - s)Y 
, 4(p + i?) 



VpTr 



[p + R)-^^^H{c{t),c{t ~s),p + R) 
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with H{c{t), c{t — s),p + R) analytic in p. This follows since — 1 = 0{z) near 
z = 0. We now substitute this back into (|3.2ip and change variables to iXp — z, 
to obtain: 



(mm 



-le 



iXR 



1 H(c(t),c{t- s),-iz/\)dz 



The integrand is analytic in A, and absolutely convergent. The power of X^/'^ 
makes the net result a ramified analytic function. The same argument can 
be applied to 71, replacing R by —R. Thus, we have shown that k'_^^{t,s) is 
analytic in A^/^. This implies that K{—iX) is analytic in A. As remarked before, 
the case K{—iX + uj) is identical, so the proof is complete. □ 

Now that it is justified, we can now write K' — Kj^. In addition, now that 
Kl{(j) and Kp{a) are defined, it is clear that Kp{a) + KL{cr) = K{a). 
Next we show that K{a) grows ad most exponentially as 5cr — > ±00. 

Proposition 3.11 K{a) vanishes as — > 00. 

Proof. We break K (a) up as K{a) = Kf{(t) + Kl{(t). The first term, Kf{(t) 
is bounded (away from cr = 0) simply by inspecting ()3.16|) . The second vanishes 
near 5tT = c» by Proposition 13.101 □ 

We have now shown that K{a) : L^{S^,dt) L'^{S^,dt) is an analytic (in 
a) family of compact operators. This allows us to construct the resolvent. 

Proposition 3.12 The operator (1 — K(a))~^ is a meromorphic (in a) family 
of bounded operators. This implies that if (1 — K{a))^^ has a pole of order n at 
a point a — (7k, we then have the following asymptotic expansion as a (Xfe.- 

(1 - Kia)r = E + Dia) (3.23) 

(ct - CTfc)J + -^ 

where D{a) is analytic nearak- Yk.j{t) solves {I — K{ak))Yk.j{t) — Ykj{t) (with 
Yk^-i(t) = 0)- The functions Ykj{t) are all L'^{SIj) functions. 

If = 0, then the same result holds, except that the poles are in the variable 
yfa instead of (cr — ijfc). 

An additional result (which we use later) is that PQy{0,t) = (1/2) J^^pg{x)dx, 
where Pq is projection onto the zero'th Fourier coefficient. 

Proof. This is merely the analytic Fredholm alternative theorem. There is 
only one technical point regarding the behavior near a = due to the fact that 
K{a) is singular there. 

This can be remedied as follows. The function y{(7,t) satisfies the following 
equation: 

(1 - Kf{ct) - KL{a))y{a,t) = y^{a,t) (3.24) 
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We expand Kp{a)y{a, t) as in the proof of Proposition l3.8l Due to the fact that 
Kp{(j) is singular only in the zero'th Fourier component (see (|3.16l) ). we find 
that: 



(1 - Kf{ct){1 - Po) - KL{a))y{a,t)+a-'/^Poy{a,t) 

= a-i/2(l/2)^Vo(x)dx + /(ai/2,t) 

Here, Pq is the projection onto the zero'th Fourier coefficient of a function. 
Take as an ansatz that Poy{0,t) — (1/2) J^ipo{x)dx. Then, since Kp{<t){1 — 
Pq) — Kl{(t) is compact and analytic in u^/^, we find that 

y(<j,t) = [1 - Kf{<j){1 - Po) - KL{<j)]-\f{<J^'\ t) 

is meromorphic in ct^/^. This implies that our ansatz was consistent. □ 



Proposition 3.13 Define K^{a) as K{a) with c{t) replaced by ec{t) (so in par- 
ticular, Ki{cr) = K{cj)). Then the position of the poles of K^{a) are ramified 
analytic functions of e, the field strength, except possibly near dk — —ioo. For 
small e, there is only one pole (Tq near the real axis (corresponding to the dressed 
bound state), and all other poles are located near a = —ioo. 

Proof. This is basically the analytic implicit function theorem, using the fact 
that Ke{a) is analytic in e, and Ko{cr) = Kp{a) (c.f. Proposition 13. 8p . 

We first show that no poles form spontaneously. Consider a compact set, 
bounded by the curve 7. Then define 

Rf^e = / [1 - K,{cr)]-~^da 

Provided [1 — K^{a)]^^ is analytic on 7, then R^y^e is analytic in e. For e = 0, 
we find that: 

-1 

/„e-™'^* (3.25) 

which has one pole on the real line, and no others. Using the fact that 
||iir(:(cr) — i<ro(o')|| < Ce (with C depending in 7), we can expand K^ia) in a 
Neumann series in e (provided 7 excludes the pole of [1 — K{a)]~'^) to show that 
[1 — K{a)]~^ is analytic for all small e. Thus = for all small 7, and is 
therefore zero when it is analytic. Since R-^^e is analytic until [1 — K^{a)]~^ is 
singular on 7, we have shown that any poles inside 7 must have gotten there by 
crossing 7. 

The same argument can be repeated to show that spontaneous poles of higher 
order do not form, except that we will need to consider R-y^e — S^fk{<^)[^ — 
K^((7)\^^da (for fkicr) a function with nonvanishing k'th derivative) instead. 



[1 - Ko{a)]- 




E 
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This further impHes that any poles which are not present for e = must 
come from tr — —ioo as e is "switched on". 

Analyticity oiak follows immediately from Theorems 1.7 and 1.8 in [221 page 
368-370] (see also the discussion following Theorem 1.7). These results show that 
any eigenvalue X{e,a) of K^{a) is analytic. Poles occur where A(e, cr) = 1. By 
the implicit function theorem, ak = <Tfc(e) is ramified analytic. □ 



Remark 3.14 If we could compute an upper bound on the norm of [1—K{a)]~^ 
as 9(7 — oo, it would be possible to show that [1 — K{<t)]^^ has only one pole: 
the analytic continuation of the bound state when E{t) = 0. The upper bound 
would make it impossible for poles to come from —ioo. In all other cases we 
have considered |i7i iSj it was possible to do this, and we believe our inability to 
do so here is a technical point rather than anything fundamental. 



3.3 Time behavior of ip{x,t) 

We have now shown that K{a) is a compact analytic operator. By the Fredholm 
alternative, (1 — K{a))~^ is a meromorphic operator family. By deforming the 
contour in (|3.6ap . we can determine the behavior of Y{t). Once this is complete, 
we can calculate $fe j(x,t) and ^M{x,t) and finish the proof of Theorem[21 

Proposition 3.15 The function Y(t) has the expansion: 

M-l rifc 

^(^) =J2Y1 akt^e"^^'Yk,j{t) + DM{t) (3.26) 

k=0 3=0 

with Yk,j(t) the residue at ak and aj.k = {2Tr/Lo){yo{ak,t)\Yk,j{t))/j\. M must 
not be greater than the number of poles of [1 — K{a)]^^. The function Dmit) 
has the asymptotic expansion: 

oo 

Dm it) ^ E D^^nt^'^^ (3-27) 

riGZ j=3 

The sum over n is convergent in P, which shows that |_Dm(^)| = 0{t^^^'^). 

Supposing that ak = for some k, Y{t) — D]\i{t) except that in p.27p the 
sum starts at j = 1 rather than J = 3. 

Proof. Because (1 — K{a))^^ is meromorphic in a, y{a,t) can be written as 



yia,t) = (1 - Kia))-^yoia,t) = ^^^.-W(^^..-(^)l^o(.,t)) ^ 

(3.28) 
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We compute Y{t) using (|3.6ap . and shifting the contour: 



Jo+ 

p-iK(M)+0+ p-iK(M)+uj 

^uj-^ e-"'*y{a,t)da + uj-'^ e'"'*y{a,t)da 

Jo+ J-iK(M) 

e"""*y(cr, t)dcr + Residues = cj"^ / e"*''*y(cr, t)da 

-iK{M)+LO- J-iK(M) 
p-iK{M)+0 

+ Lo-^ / e-'"*2/(a + 0+, t) - e-^'^+^-^'y{(j + LO-,t)d(T 

Jo 

+ Residues (3.29) 

If [1 — K{a)]^^ has more than M poles, then we make K{M) sufficiently large 
to collect M of them; otherwise, we simply collect all the poles. The residue 
term is given by: 

Af-l rifc 
k=0 j=0 

stemming from the M poles with 3(Tfc > —K{M). By p.6cp . we can change the 
integral in the second to last line of (|3.29p to: 

l—iK(M) 

uj-^ e^'^^yia + 0+,t) ~y{a + 0^,t))d(T (3.30) 

Ja 

Note that y{cr,t) is analytic in a^^^, and thus y{a + 0+,t) — y{<7 + 0_,t) can 
be expanded in a Puiseux series in cr^/^ (and a Fourier series in t). Watson's 
lemma yields: 



p — ioo 



neZ j=0 

oo 

~ CO-' J2 e""' E D,,nT{j/2)t-^/^ (3.31) 

nSZ j=3 

This is what we wanted to show. 

When (Tfc = 0, the result follows simply by noting that the sum over j in 
p.30p starts from j = —1 rather than j = 0, thereby letting the sum on the 
right of (|3.30p start a.t j = 1 instead of j = 3. 

The integral from —iK{M) to ~-iK{M) + uj decays at least as fast as 
0(e-^(*^)*), and is included in Dk{t). □ 

We now reconstruct ip{x, t) in the velocity gauge. The basic idea is as follows. 
We know that ipy{c{t),t) = Duit) + e-''^>'*Ykj{t). Using the fact that S{x ~ 
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c{t))ipv{x, t) = S{x — c{t))ipy{c{t), t), we find that ^l^v{x, t) satisfies the following 
equation: 

idtipv{x, t) — —d1i]jy{x, t) — 25{x — c{t))tpy{x, t) 

= -dliP^ix.t) - 25{x - c(t))^^(c(i),i) 

= -dlM^, t) - 25{x - c{t))[DM{t) + e-**'=*rfc,,(t)] 

Proposition 3.16 The expansion (|1.7p holds. 

We first state a result, proved in Appendix [B] which is necessary for the 
proof. 

Proposition 3.17 Let G{a) be the Green's function for the equation: 

{a + idt + dl - b{t)d^)u{x, t) = -2S{x)f{t) (3.32) 

so that for Scr > 0, u{x,t) decays as x ^ ±og. Here, u{x,t) is in L^(55) for 
each x, and fit) e L'^(5^). Then G{a) can be analytically continued to the 
region 3cr < 0. The function u{x,t) — G{a)[—2S{x) f (t)] has the expansion: 

Am,± — tWo' + riiuj (3.33b) 

where f{t) i— » {um,R,Um.L\ is a mapping from L'^{S^) — > x {L,R}). G{a) 
is also a continuous map, analytic in cr^/^ from L'^{S^) — > LP'{Br x S}^) with 
Br ~ {x : \x\ < R} for any R. Near a = 0, the we have G{cr)S{x)f(t) = 
(j-i/2(^2/2)Po/(i) + 0{1), with Pof{t) the projection onto the zero'th Fourier 
coefficient of f(t) and the 0{1) term being analytic in a^^^. 

Proof of Proposition . 1 61 We work in the magnetic gauge, to simplify this part 
of the problem. Note that ipBix^t) = ^pv{x + c{t),t), so in particular, ^pB{0,t) = 
ijjv{c{t),t) — Y{t). Moreover, recall that the Zak transform commutes with 
periodic operators, such as the coordinate transform (x,t) i— > (x + c{t),t). 

Additionally, in what follows, the notation 0(cr^/^) denotes a function ana- 
lytic in (T^/^ taking values in L'^{S^, dt). 

By Zak transforming the Schrodinger equation in the magnetic gauge, we 
obtain the following (with ^'((t, a;,i) = Z[ili\{a,x,t)): 

{a + iat)*(cr, X, t) - Vo(a;) = [-A + 2ih{t)d^] ^{a, x, t) - 2S{x)^{a, 0, t) 

= [-A + 2i6(t)5J ^{a, X, t) - 25{x)y{a, t) 

Bringing all terms besides —25{x)y{a,t) to the left, the initial condition to the 
right and inverting the differential operator yields: 

^{a, X, t) = [+a + idt + A- b{t)d^r^ M^) 

- [+cr + idt + A^ b{t)d^]''^ '2S{x)y{a, t) (3.34) 
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The second term is given by —G{a)2S{x)y{(T,t). To compute tp{x,t), we need 
to invert the Zak transform. Thus: 



Jo 

JO J-iM J-iK{M)+LU 

+ Residues (3.35) 

Note that the first term of (|3.34p can be equivalently written as 

Z[{a + idt + dl)-^Mx)\{x + c(i), t) = Z[e'^'*Vo(a:)](a, x + c{t), t) 

By Proposition [Sill for each x, this term takes the form (l/2)cr^^/^ / il'Q{x)dx + 
f{a^^^,t) with f{a^^^,t) varying with x. Thus, the first term is bounded by 
(27gC|Q<T| fgge Proposition l3.6p . The second is given by —2G'((t)(5(x) j/((t, t). The 
integral along [—iMj—iM + uj] decays like 0{e~^*). Thus, (|3.35p becomes: 

ip{x,t) 

l—iK(M) i—iK(M)+uj 

= I e-"'^^{a,x,t)da - I e"""*^'(cr, x, i)dCT + Residues + 0(e"*^*) 

Jo Jul 

(3.36) 

We show that the contour integral in p.36p gives rise to the dispersive part, 
while residues give rise to the resonance. 
The Dispersive Part 

To compute the integral term of (|3.36p . note that we must compute: 
*M(a;,<)=/ e"*'"**(cr,x,t)do-- / e''"'*'i>{a, x,t)da = 

Jo Jul 

e-^'"*Z[e*^="*Vo(a;)](a, a;+c(i), i)dcr-e~'('^+")*Z[e*^=^Vo(a;)](CT+t^, a;+c(t), i)da 

p-iK{M) 

~2 e'"'*G{a)6{x)y{a, t) - e-'^''+'^^' G{a + uj)5{x)y{a + t)da (3.37) 

Jo 

Since Z[f]{a + uj,t) = e*"*Z[/](cr, i), we find that: 

e-''"*Z[e''^'Vo(a;)] (a, x + c{t), t) - e-*('^+'^)*Z[e*^' Vo(a;)](CT +io,x + c{t), t) 

= e-*'"*Z[e'^' Vo(a;)] (cr, x + c{t), t) - e-*'^*Z[e^^' Vo(a;)](a -0,2; + c{t), t) 

(3.38) 
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Using the fact that 2:[e*^'Vo(a;)](<5, a;+c(0, i) = fT-i/2(l/2) / iljo{x)dx+0{(j^/'^) 
(by Proposition 13. 6p we find that: 



Plugging this into p.37p yields: 



(15371) = i^^''^ j Mx)dx- 



-iK{M) 







l.-iK{M) 

-2 e-"'*G{a)6{x)y{<T,t)-e-''^''+'^'>'G{<T + uj)5{x)y{(j + uj,t)d<T (3.39) 

Jo 

Again using the identity Z[f] (cr + u!,t) = e*'^*Z[/] {a, t), we find: 

-iK{M) 

lO 

l—iK(M) 

-2/ e^*''*[G(f7 + 0)(5(x)y(cr + 0,t)-G(a-0)J(a;)2/((7-0,t)]dcr (3.40) 
Jo 

Since G(cr)^(x) = a-i/2(l/2)Fo + 0(1) near cr = 0, we find the second integral 
term in p.40p becomes: 



= t'^'^ J Mx)dx + J e~'^'Oia^^^)da 



-iK(M) 



{a + 0)-'/\l/2)Poy{a + 0,t) + 0{l)y{a + 0, t) 

da 

a-^^^Poy{0, t) + e-'"*0(ai/2)da (3.41) 



(a - Q)-^'\l/2)PQy{a, t) + 0{l)5{x)y{a - 0, t) 

r-iK{M) 



Using the fact that Poy{0,t) — (1/2) J il;o{x)dx (see Proposition I3.12| ). and 
plugging into yields: 



'fM{x,t)^t-^/^ I Mx)dx + J e-"''0{a^''^)da 



'iK(M) 

2 / e-*'"*(l/2) 
'o 



'ijji){x)dx 



da 



-iK(M) 



e-'"''0{a^^^)do 



l—iK(M) 

/ e-'"*0(ai/2)da 
Jo 
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This is the Laplace transform of a function which is analytic in cr^/^, which by 
Watson's lemma yields (ll.9[) . 
The Residue Term, o-fc 7^ 

By substituting p.23p into p.35p . we find that when 7^ 0, the residue 
term (for each pole) takes the form: 

with gfc^j ^ (yfcj(t)|j/o(o-, t)) and <^k,j{x,t) = G{dk) 5{x)Y k^j{t). Thus, by Propo- 
sition [3II3 (in particular (|3.33ap ). we have proved (jl.lip . 

This imphes that \1/((t, cc, t) has a pole &i a = dk with residue <^k,j{x, t). Since 
\l/(cr, 0, t) = y{a, t), and y{a, t) has a pole at cr = d-fc with residue Ykj (0)t, we find 
(equating the poles) that $fcj(0,i) = >fej (t)- Thus, (5(a;)yfej(t) = (5(a;)$fcj (0, t), 
which implies that: 

$fcj (a;, t) = -2[(Tfe + + A - 5(i)a:.]^'<5(2;)$fcfej (0, t) 

Applying [dk + i9t + A — b{t)dx] to both sides and rearranging yields (jl.Sp . 
The Residue Term, o-^O 

Supposing CTfc = (for some k), by Proposition 13. 131 we find that y{<T,t) = 
Ykj{t){Yk^.j(t)\yo{0,t))a-^/'^ + D{a,t)yo{a,t). By PropositionEH we find that 
Poyia,t) = (1/2) J^'tljo{x)dx; thus we find that PoYkj(t) = 0. 

Since G{a) — a~^^^Po + 0{1), and Poy{0,t) = 0, this means that for small 

a: 



The "residue" term therefore becomes: 



a- / e-'"'G{G)6{x)Yk,,{t)^ 
w Jo V 



The 0(t~^/^) term comes from computing a Laplace-like integral of a Puisseux 
function, and can be incorporated into ^M{x,t). The 0{t~^^^) term is the 
zero-energy resonance. This completes the proof. □ 

We have thus far proved all of Theorem [2] except for (jl.lOp . 

Proposition 3.18 Suppose that 3(7^ — 0. Then il.lO\) holds and ij^n '^ — 

all n < 0. Furthermore, the pole is of order 1. 
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Proof. It is clear that unitary evolution implies: 



^ , , \ip{x,t)\ dxdt<l (3.43) 
27r Jo J-R 

If the pole is of order greater than 1, then: 



E 



But the second two terms decay, while the first grows with time. This contradicts 
unitary evolution, unless Uk = 0. Thus the pole must be of first order. 

Now suppose that in the expansion of ^kfi{x, t), at least one V'n'^ with 
n < 0. Then ^k,o{x,t) will oscillate with x rather than decay. This implies 
that: 



2^ X 1-. 



\^k,j(.x,t)\^ dxdt > CR 



-R 

for sufficiently large R and some C > 0. On the other hand, the rest of %l}{x^t) 
(the dispersive part, and the exponentially decaying poles) which we denote 
R{x,t) decays with time. This implies that for t > tn (with Tr chosen large 
enough so that \R{x,t)\ < e/y/2R that: 

mx,t)\\ = \\^kfi{x,t) + R{x,t)\\ > \\^kfi{x,t)\\ - \\R{x,t)\\ > VCR-e 

Selecting R> {2 + e)/C causes ||'!/'(x, f)|| > 1, contradicting unitary evolution. 

Intuitively, what this means is the following. The modes V'n'^ with n < 
correspond to radiation modes. If such a mode is nonzero, then (^kfiix^t) will 
be emitting "radiation" without decaying, which is clearly impossible. □ 



4 Concluding Remarks 

In this paper we studied the interaction of a simple model atom with a dipole 
radiation field of arbitrary strength. We obtained a resonance expansion, in 
which resonances can be resolved regardless of their complex quasi-encrgy. In 
particular, we obtained a rigorous definition of the ionization rate 7 = —2'^dk 
and Stark-shifted energy, Jfd-fe for the fc-th resonance. We further showed that 
complete ionization occurs (7 > 0) when E{t) is a trigonometric polynomial. 
Some possible future directions of research include: 

4.1 Perturbative and numerical calculations 

The main feature of our method is that it turns a time dependent problem 
on M into a compact analytic Fredholm integral equation. This implies that a 
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family of finite dimensional approximations can be used (in the Zak domain) to 
approximate solutions to the time dependent Schrodinger equation. 

We believe that the quasi-energy methodology used here and in related pa- 
pers [51 [HI Uj can be used for quantitative calculations of realistic physical sys- 
tems. Perturbative calculations along these lines have recovered Fermi's Golden 
Rule and the multiphoton effect. 

4.2 Resonance theory 

Significant effort has been devoted to the rigorous definition of resonances and 
quasimodes. The best results we are aware of are those of [THl [31] , based on 
complex scaling, and those based on analytic continuation of the S-matrix, e.g. 
[T]. We provide an alternative definition: a quasi-bound state is the coefficient 
of an exponentially decaying term in the asymptotic expansion for 'ijj{x, t) near 
t — oo. We hope to use this definition to provide a more complete picture of 
the time evolution of ip{x^t). 

4.3 Extension to 3 dimensions 

In the case of ifo = —A — 2^(2?") with x G M^, a similar equation to (|3.4p can be 
derived. Due to the fact that 5{x) is not in H~^{E?)^ ip{x,t) becomes singular 
at t = 0^ . This can be remedied by considering weak solutions, and an equation 
similar in most respects to (|3.4I) can be derived which governs the evolution [T3] . 
For this reason, we believe most results can be adapted to this case, as has been 
done for Ha = -A^ 2d{x) + E(t)S{x) [51[5]. 

A Proof of Proposition 12.11 

We observe that by the results of Section [3l if a bound state exists, then: 



Setting z = e and y{z) — Yk{t), we wish to show that y{z) = f{z) + g{z) 
with /, g both entire of exponential order 2n. This is equivalent to showing that: 



\Yk{t + ia)\ < Ccxp[C"exp(|2iVwa|) 



The function Yk{t) satisfies the equation: 




with k'(t,s) as defined in (|3.19bp . Thus we obtain the bound: 




\k' (t + ia, t + ia — s) \ \Yk{s)\ ds 



(A.l) 
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and it suffices to bound \k'{t + ia,t + ia — s)\. Prom the definition of k'{t, s), 
we find: 



k'{t + ia,t + ia — s) 
= -[ - 



, t:,ujp-\-a—iuj(t— s) 
R+Oi J- — e ^ ' 



exp 



{c{t + ia) - c{s)f 
Ap 



dp 

VP 



Supposing a/uj > 1 (we are interested in the behavior as a ^ oo), the integrand 
is analytic for z — re'^, < r < 1 and Q < 6 < tt. Thus, we can deform the 
contour from M + Oi to 7 = d{z : 9z < Oor < 1}. 

Note that for some constant C, \c{t + ia)\ < Ce^'^^"^, since c{t) is a trigono- 
metric polynomial of order N . 

We find that there are three regions of integration which contribute to k'{t + 
ia,t + ia — s). The regions of integration contributing come from the region 
near 1 — e'^'P+«-»'^(*-«) = Q (the pole of the integrand), large p and small p. 

If the pole is closer to M than tt/uj, we deform 7 up to encircle it at a distance 
piLo. Otherwise, we ignore it. Therefore, in any case, for 2; € 7, l — e'^P+"~*'^(*~*) 
is uniformly bounded away from zero. 

We then split 7 = 7< U 7> U 7^ where 7< = {p e 7 : |p| < (Ce^"l"l + 
lk(s)||j;,oo)^} and 7> = 7 \ 7<. We therefore find that: 



\k'{t + ia,t + ia — s)\ < |residue| 



C 

+ c 



2 gcjp+a— iaj(t— 5) 



exp 



/ 

"'7> 



gOJp+a— ia;(t— s) 



exp 



{c{t + ia) - c(s))2 
4p 

{c{t + ia) - c(s))2 
Ap 



The residue can be bounded by: 

I residue I 



< C 



,<7{—a+ko{t—s))/ui 



exp 



(c(i + ia) - c(s))^ 



dp 



dp 



< C 



4{-a + iu{t- s))/u;J J V{-a + iuj{t - s))/uj 
< C exp{C \c{t + ia)^^) < Cexp(Cexp(2iVu; |a|)) 



We bound the integral over the compact region 7< simply by taking absolute 
values: 



^ gOJp+a— iaj(t- 



exp 



(c(i + ia) - c(s))' 



4p 



dp 



< |7<|Cexp(Cexp(2A/'w|Q;|)) 
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For the integral over 7>, we use the fact that if|z|<l, |e^~l|<e|z| 



7> 



exp 



{c(t + ia) - c(s)f 
(Ce^-I"l + |jc(s)|j 



dp 



^ „ QUjp-i-a~iLo{t — s) 



\p\ 



dp 



2 gojp+a— ic(j(i— s) 



P 



-3/2 



dp < Cexp(Cexp(27Vw |a|)) 



Combining these estimates, we find that k'(t + ia, t + ia — s) has the required 
growth as a ^ oo, hence Yfc(t) does. The same argument apphes as a ^ —oo. 



B Proof of Proposition 13.171 

We state a few results we need. 

Theorem 4 (T. Kato, [22, page 368]) If a family T (a) of closed operators 
on X depending on a holomorphically has a spectrum consisting of two separated 
parts, the subspaces of X corresponding to the separated parts also depend on a 
holomorphically. 

Remark B.l A few words of explanation are in order. In [52], they are given 
in the commentary following the theorem. 

The analytic dependence of the separated parts of the spectrum means 
the following. Let M^r, M'^ be the spectral subspaces of T{a), related to the 
two separated parts. Then there exists an analytic function U{a) (called the 
transformation function), with analytic inverse, so that — U{a)MQ and 
M'^ = U{a)Mlj. For fixed a, both U{(j) and U~^{a) are bounded operators on 
the Hilbert space. 

In addition, the spectral projections Pm{o') and Pm'^ct) can be written as: 

PMicr) = U{a)PMiO)U-\a) (B.la) 
Pm{(j) = U{a)PM'{0)U-\a) (B.lb) 

We now prove a Lemma which allows us to reconstruct ^'(ct, x, t) given solely 
information about ^'(ct, 0, t). The basic idea is to treat the Schrodinger equation 
as an evolution equation in x, with a "Hamiltonian" that is periodic in t. 

Lemma B.2 Define the Hilbert space H = A{Sl,dt)®L'^{Sl,dt), with A{Sl,dt) 
defined by the norm: 

\\I{t)\\\^{fm^ + ^dt\f{t))L^siM) (B.2) 

Note that the norm on ^(5^, dt) is equivalent to the norm on H^/'^{Sl,dt) except 
at a = 0. Then there exists a sequence with < inf„j |A'^„j| < sup„j \Nm\ so 
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that if we define 4>m,± o,nd \m.± os follows, 

4>m,± = Nr 



9-1/2 \-l -imwt =FA,„,±c(t) 



Am,± — Ti^/cr + mijj 
then (j3m,± "is a Riesz basis for TL. Furthermore: 



(B.3a) 
(B.3b) 



H 



1 

a + idt 





b{t) 



is diagonal in this basis, with H(j)m.± = ^m.±4'm,±- 

Moreover, if we define Ti.'^ as the span of {4>m,+}mei o,s the span 

of {cj>rn,-}m£Z then e^^ is defined, bounded and analytic in a on Ti^ for x <Q, 
and on 'H~ for a; > 0. 

Remark B.3 Note that the norm on A{Slj,dt) controls the Xq\. term in the 
first component oi (t>Q,±, which would otherwise blow up as ct ^ 0. 



We are nearly ready to prove Lemma IB. 41 First a minor technical point 



Remark B.4 Consider the sequence y/a + nuj, with 3?cr S {0,lu). For n nega- 
tive, 3^/(7 + nuj grows like ^/\rl\■ For n positive, a + nuj — 0(n^^/^), and is 
uniformly bounded below. 



Proof of Lemma \B.2\ It is a simple calculation to show (f>m,± are eigenvectors 
of H with eigenvalues Am.±. To show that is a Riesz basis for "H, we show 

that iJ is a bounded perturbation of a normal operator. Consider the family of 
operators (analytic in C) on Ti.: 



u 



1 

a + idt 



u 





Cb{t) 



u 

Urn 



Consider also the family of vectors (parameterized by Q : 



0-1/2 \-l -imujt C^m.±c{t) 



Njnx is a normalizing constant which is defined implicitly; we discuss it 
below. For C = 0, N„,x = 1- 

A simple calculation shows that (0^^,Am,±) are eigenvector /eigenvalue 
pairs of Hq. In particular, each Xm,± is separate from all the others. For 
C = 0, they are also orthonormal in Ti.. Let P,^(C) be the associated spectral 
projection operators, given by 



(ijc - zy^dz 



32 



where 7,„ is a closed curve containing only ^, and no other eigenvalue of iJ^. 

Let C/(C) be the transformation function of Theorem [4] (on pagelSTl see also 
Remark IB. II and (|B.ip V Since each eigenvalue Xm,± is separated from all the 
others and varies analytically (except near ct = 0), Theorem 2] implies that: 

Note now that P^(0) — ( • \4>m.o)'f'm.o- Therefore, we can write: 
We know that U {C)4'm o a vector in the direction 

2-l/2g-imwtg=FC-^m,±c(i) 

but this determines U{C)4>m o only up to a constant. Nrn,.(; denotes this constant 
(see (|B.ip in Remark lB.ip . Because U{Q is bounded above and below, we find 

that that < ||[/(C)"^|p^ < W^xl < 11^^(011 • Note also that N^x is not 
necessarily real. 

To compute the expansion of a function ?/;(<) in this basis, we use the formula 
tp^ = {U{l)~^ip\(j)'^ o). Since (j)^ q is an orthonormal basis, this set of coefficients 

is clearly in with P norm bounded below by ||t/(l)||~^ 11-0(0, and above 
by||C/(l)-i||||V(0,i)||^. 

Finally, we need to show that J2m^miO + -Pm(C) — 1; where the sum is 
interpreted to converge in the strong topology. The sum is strongly convergent 
when C = 0; since the ^^(0) are orthogonal projections. Now multiply on the 
left and right by U{Q and U{Q~^ (recalling Remark [BT] and (jB.ip ). which are 
continuous operators (in norm and therefore in the strong topology): 

This completes the proof of the Riesz basis property. To show boundedness of 
e^^, simply note that the real part of the eigenvalues of H is bounded above 
on and bounded below on 'H~ (though not uniformly in a). Thus, e^^ is 
bounded on 7i+. Analyticity follows simply by observing that the eigenvalues 
and eigenfunctions are analytic in cr^/^. □ 

We are now prepared to prove Proposition 13 . 1 71 
Proof of Proposition \37T7\ Note that p.32p can be rewritten as: 

dx 





)-"{ 




V ux , 
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Away from x = 0, the solution u{x, t) can be written (formally) as: 

At a; = 0, the two matching conditions need be satisfied: 

m(0+,<) - u(0",t) = (Continuity) 
dj,u{id+ ,t) - d:ru{0^ ,t) = -2f{t) (Differentiability) 

For Q(T > 0, Xm,+ always has positive real part and A„i._ always has negative 
real part (recall (jB.3bp ). Thus, if u{x, t) is to vanish as a; — > ztoo, we find that: 

u{0-,t) 



u(0+,t) 



Um,R4>m-{t) 

m 

y^Mm,L0m,+ (O 
m 

Since 4>m,± is a Riesz basis and [0, f(t)] G 7i, we can write: 

fit) j = ^f"-'^,+ ^''n,+ + fm,~-(t)rn- (B.5) 
^ m 

Choosing Wm.i? = fm.- and Um^^ = ^fm.+ solves p.32p . at least on a formal 
level. Since w(0+,t) G and M(0~,t) S , (|B.4p makes sense. Since e^-'^ is 
bounded and analytic provided |a:| < i?, this is thus an analytic mapping from 
L^iSl) ^ L^iBn X Si). 

Now observe that both (|B.4p and (jB.Sp can be analytically continued in 
(T, and the continuation also solves ()3.32p . therefore G{(j) can be analytically 
continued in a as well. 

We now need only determine the behavior near a — 0. By Taylor-expanding 
(|B3a)) in cr^/^ we find that: 

\ / o(i) 



while 



Thus, near cr = 0, we find to leading order (plugging (|B.6p into (jB.sp ) that: 

/m,+ fm,— 

A^o2-i/'(/,„,+ + /,„,-) = Po/(i) 

with Pof{t) projection onto the zcro'th Fourier coefficient. This implies that 
fo,± = uoM ^ -uox = 2i/2iV-i(l/2)Po/(t)+0((Ti/2). On all other coefficients. 
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the behavior is analytic in a since A„i,± is analytic in a for m 7^ 0. Thus for 
small a: 



and therefore G{a)5{x)f{t) = (l/2)cr"i/2[Po/(t)] + 0(1) near cr = 0. □ 



C Wellposedness 

Given that Y{t) exists and is smooth (easily seen by using a Banach fixed point 
argument on (|3.4p ). we need to extend Y{t) to ipv{x,t). This is done by means 
of p.2p : the main thing to show the extension is in L^(R) for each t. Since 
ipy^o{x, t) is in L^, we need only show that: 

„ /■* f iix - cit ^ s)f\ , , , , ds 
2i / exp ^ ] '-^ i;,{c{t -s),t- s) 



\ 4s y y/Anis 

= 2i f e"t''^'''''''^Y{t~ s)ds (C.l) 
Jo 

is in L'^{W,dx). This is done by stationary phase. The phase of the integral 
(incorporating the s^^^^ from the integrand) is (f>{x,s,t) = (x — c{t — s))^/s + 
i(l/2)lns. The phase is stationary when: 

f>M A H^-c{t-s)f , -2{x-c{t-s))b{t-s) I 

Provided x > sup^ \c{t — s)|, we find that ds<j){x, s, t) > 0, and thus the phase 
is never stationary. Further, note that ds4'{x, s,t) — 0{x^) near x — 00. The 
integral over [e,t] becomes: 

2i f e'''^("'"'*)r(t - s)ds = 2 / i9,0(x,s,t)e^'^("^"'*) ( ^-77— ^] ds 



= 2 



e 



....t) ( Yjt-s) 



s=t 



^ds(t>ix,s,t) ^ 

_ 2 r ^r4>[...s,t) d.^i^, t)dsY{t -s)- Y{t - s)dl4>{x, s, t) 

Je {ds<j){x,S,t)f 

Since e^'^^^'^'*^ — 0(s^^/^) near s = while ds<l){x,s,t) = 0{s~^) near a; = 0, 
we find that e*'^^^''*'*-'/^s0(2:, s, i) ^ as s ^ 0. This implies that the first term 
on the right of (|C.2[) is bounded even as e ^ 0. A simple calculation shows that 
dg(j){x, s, t)/{ds4>{x, s, i))^ behaves like s; plugging this into the integral term of 
(jC.2p shows the integrand behaves like s^/^ near 5 = 0. Since the integration 
is of a bounded function over a compact region, and the bounded function is 
0{x~'^), we find the integral decays like x~'^ as well (even as e — > 0). 
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The same argument can be applied to compute the x or ^-derivative of this; 
the main difference is that we must replace Y{t — s) by dx4>{x, s, t)Y{t — s) or 
dt(t>{x,s,t)Y{t — s), which causes the integrand in (|C.2p to behave like s~^/^ 
near s = 0. This is still integrable0. 
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Figure 1: A plot of K£(2;), for the specific choice of (t{z) = — 5z—5z~^ + z~'^. 
As is apparent from the figure, the regions ^€{z) > and ^<^{z) < are 
approximately sectors for sufficiently large \z\. 
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